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SUMMARY 


Estimates  of  critical  distances,  up  to  which  dust,  aerosols,  and 
(toxic)  fumes  may  be  driven,  under  the  influence  of  various  wind  and 
diffusivity  conditions,  require  that  three-dimensional  solutions  to  the 
problem  of  wind-driven  contamination  be  derived,  for  application  to 
range  safety  problems.  The  so  called  "conservative,"  two-dimensional 
solutions  (which  might  be  of  seme  value  in  certain  meteorological 
situations),  are  in  general,  by  orders  of  magnitude,  too  large  to  be  of 
practical  value  (see  part  III,  page  5).  After  all,  nature  is  three 
dimensional;  this  feature  determines  the  rate  at  which  the  concentration 
falls  off  in  the  downwind  direction. 

In  part  III,  analytic  solutions  for  constant  wind  and  constant 
diffusivity  have  been  derived.  By  using  the  method  of  separation  of 
variables,  the  three-dimensional  concentration  functions  resulting  from 
two-dimensional  as  wall  as  three-dimensional  sources  are  deduced.  In 
either  case,  the  intensity  q(x,  0)  falls  off  with  the  first  power  of  the 
downwind  distance  x,  as  expected;  but  the  strength  factor,  which  depends 
on  the  dimensions  of  the  source  in  the  x-,  y-,  and  z-directions,  offers 
some  new  insight  (3ee  figure  l).  Even  in  three  dimensions,  the  resulting 
critical  distances  are  too  large  to  be  useful,  so  long  as  the  assumption 
of  constant  U  acd  K  is  retained  (see  part  III,  page  5). 

In  part  IV,  three-dimensional  analytic  solutions  have  been  derived, 
under  the  assumption  that  wind  speed  U  and  eddy  diffusivity  K  vary, 
either  in  accordance  with  the  conjugate  power  laws,  or  in  a  more  general 
fashion.  The  results,  following  therefrom,  can  be  made  to  agree  with 
published  observational  data;  in  particular  with  the  general  observation 
that  Che  rate  of  decline  of  concentration  with  the  downwind  distance  x 
and  crosswifid  distance  y  along  the  ground  surface  obey  a  certain  power 
las-/. 


In  section  one  of  part  IV,  a  new  method  of  solution  is  suggested 
which  may  be  used  for  solving  diffusion  problems  of  a  more  general 
nature:  a  "method  of  superposition  of  parabolic  sources,  sinks,  »tc.u 

In  the  remainder,  five  particular,  three-dimensional  parabolic 
source  solutions  have  been  derived. 


The  five  mentioned 

solutions 

are  characterised  by: 
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and  Kg  are  the  diffusivities  in  the  y-  and  z -direct ions,  p  is  the 
exponent,  which  indicates  the  rate  of  decline  of  concentration,  with 
downwind  distance  x,  The  examples  are  ordered  after  Increasing  p„ 


PART  III 


THE  CASE  OF  CONSTANT  WIND  VELOCITY 
AND  CONSTANT  DIFFUSIVITY 

THREE-DIMENSIONAL  DIFFUSION  FOR  A  TWO-DIMENSIONAL  SOURCE 

In  part  I,  the  writer  has  derived  an  analytic  solution  of  the 
problem  of  diffusion  of  wind-driven  _cr.taminat.ion,  in  two  dimensions; 
wind  velocity  and  dlffusivity  were  assumed  to  be  constant.  In  part  III, 
an  analytic  solution  of  the  analogous  problem  in  three  dimensions  will 
be  derived. 

As  in  figure  1-1  of  part  I,  the  wind  velocity  U  is  in  the  direction 
x,  and  the  source  is  located  in  the  y-z  plane  (x  “  0)  and  may  extend 
between  0  <  z<  i ,  -  ra  <  y  <  J-  m.  (Later  on,  we  shall  assume  a  three- 

dimensional  source.)  We  desire  the  three-dimensional  distribution  of 

the  concentration,  q(x;  y,  z) ,  which  has  to  satisfy  the  differential 
equation: 

M  .  K  is  +  all  u) 

3x  U  Qy^ 

Again,  we  stipulate  the  boundary  conditions: 


—  A  jk 

q  °  qc  (y,  z)  =1  for  -l  <  z  <  +  and  -  m  <  y  <  +  ra 

(2) 

=  0  outside  of  this  rectangle  in  the  y-z  plane 


The  symmetrical  form  of  equation  (2),  about  the  x-y  plane,  guarantees 
that  the  second  boundary  condition; 

K  2a  -  o,  for  z  *  0,  any  x,  any  y  (3) 

oz 

will  be  fulfilled,  which  signifies  the  impenetrability  of  the  ground 
surface. 
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Figure  1.  Location  cf  Two-Dimensional  Source. 


We  assume  a  solution 


q(x;  y,  z)  =  q1  (x,  y)  •  q2  (x,  z) 


then 
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Thus,  the  form  of  and  q2  is  as  in  two  dimensions;  as  a  consequence, 
not  only  differential  equations  (1),  but  also  boundary  conditions  (2)’ 
and  (3)  are  fulfilled. 
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The  solution  of  the  problem  can  be  written  in  the  form 


+/  +m 


Evaluating  either  factor  of  equation  (6)  as  in  part  I,  we  get 

q(:c;  y,  2)  =  I  +  J  (u^  +  I  (u2)  *  i  +  I  (wx)  +  I  (w2)  (7) 
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Here 


By  using  the  first  term  in  the  asymptotic  expansion  of  the  erro 
integral,  we  may  write;  in  the  vicinity  of  che  x-axis,  and  for 


In  three  disens ions 


and  for  &  "  ta: 
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The  concentration  at  a  given  distance  x,  is  proportional  to  the  area  of 
the  source  and  the  strength  of  the  wind  but  inverse  to  the  diffusivity  K, 
which  diffuses  the  contamination  in  the  y-  and  z-directions . 

APPLICATION  TO  RANGE  SAFETY  CONDITIONS 


Let  us  assume  that  the  critical  contamination  value  q  is  considered 
to  be  1  per  cent  of  the  ’'initial"  contamination  qQ  =  1,  at  the  source. 

We  would  like  to  know,  at.  what  distance,  xcr,  the  value  qcr  is  reached, 
particularly  downwind  on  the  ground.  Depending  on,  whether  we  use  the 
two-  or  three-dimensional  model,  we  get: 
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Strip  source: 
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(b)  Area  source: 
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Th'?  data  given 


In  contrast  to  this  two-dimensional  result,  we  compare  the  result, 
obtained  in  three  dimensions;  equalizing  equations  (11)  and  (11a) 
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If,  according  to  the  above,  we  assume  i.  »  m  =  30,  ~  =  15  m“^ ,  we  get: 


(u} 


1.2  •  10: 


The  three-dimensional  distance  is  obtained,  by  taking  the  square  root  of 
the  two-dimensional  distance,  ar.d  in  addition,  multiplying  with  the 
strength  width  in  the  y-direction  (ratio  of  wind  width  and  dif fusivity) . 

THREE-DIMENSIONAL  DIFFUSION  FOR  A  THREE-DIMENSIONAL  SOURCE 


If  the  source  has  also  an  extension  in  the  x-direction,  -  d<x<+  d, 
(see  figure  2)  then,  if  we  skip,  for  the  present,  the  y-direction  we  get 
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Set  x  +  0  •  x,  then 
1  x+d 
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x-d 


Expanding  the  square  roots  into  polynomial  series,  we  have 
q~2  d  -  (0  •  k-  jl  +  I  + - ( 

4  W  8  x2  j 

To  get  consistent  results,  we  have  to  use  also  the  asymptotic 
expansion  of  the  error  integral  up  to  the  second  term: 
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Now,  evaluating  the  integral,  and  observing  that 
1  1  1 
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we  obtain  finally: 


(18) 


The  first  term  of  the  series  indicates,  that  the  source  has  gained 
strength  by  a  factor  2d.  The  second  term  takes  into  account  that  the 
initial  contamination  is  distributed  in  the  x-direction  (multipole  effect) . 

Analogous  results  are  obtained,  if  we  add  to  equation  (18)  the 
results  obtained  for  the  y-direction: 


PART  IV 


WIND  SPEED  AND  DIFFUSTVITY  VARYING  ACCORDING 
TO  THE  CONJUGATE  POWER  LAWS 

INTRODUCTION 

In  the  following,  the  idealization  that  wind  speed  U  and  diffusivity 
K  are  constant,  will  be  discarded.  Since  it  is  well  known  that  at  the 
earth's  (or  ocean)  surface  the  wind  speed  must  decrease  to  zero, 
accord ing  to  the  physical  laws  of  viscuous  flow,  and,  since  the  air  in 
the  atmosphere  is  in  turbulent  motion,  a  natural  assumption  is  that  the 
variation  of  wind  speed  with  altitude  obey  the  y  power  law  of  turbulent 
flow,  as  derived  by  von  Karman  and  confirmed  in  numerous  laboratory 
tests.  At  high  altitudes,  wind  conditions  will  vary  according  to 
orographic  conditions,  and  finally  there  may  be  an  inversion  layer. 

Such  complicated  conditions  can  be  best  dealt  with  by  numerical  integra¬ 
tion  of  the  differential  equation,  for  which  an  example  is  given  in 
part  II.  Representative,  simple  conditions  should  be  treated  analyti¬ 
cally,  when  possible.  This  is  done  in  the  subsequent  part  IV. 


In  marked  contrast  to  the  case  of  constant  U  and  K,  where  we 
obtained  a  range  of  concentration  characterized  by 

1 
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q(x,  0) 


in  two,  and 


q(x,  0)»=(^Jin  three  dimensions,  we 


shall  now  expect  a  decrease  in  q(x,  0)  with  higher  negative  powers  of  x. 
In  the  case  of  constant  U,  the  wind  velocity  at  the  ground  is  larger  than 
zero,  and  hence,  the  emanating  contamination  is  carried  farther.  For  a 
speed  distribution  following  the  y  power  law,  the  ground  layer  is  dragged 
by  turbulent  mixing  with  the  upper  layers,  which  are  moving  with 
increasingly  greater  speeds.  Therefore,  the  highest  concentration  is 
close  to  the  source  but  falls  off  more  rapidly  with  distance  downwind 
than  the  concentration  for  the  constant  wind. 

SUGGESTION  OF  A  NEW  METHOD  OF  SOLUTION 


In  part  III,  the  writer  has  derived  three-dimensional  solutions  to 
the  problem  of  wind-driven  contamination,  whereby  it  was  assumed  that 
wind  speed  and  diffusivity  are  constants.  "The  failure  of  these  solu¬ 
tions  tc  conform  to  the  meteorological  observations,  produces  unassail¬ 
able  evidence  that  eddy  diffusion  in  the  atmosphere  cannot  be  represented 
by  the  'Fickian'  equation,"  v.i.e.,  the  diffusion  equation  with  constant 
U  and  K) .  Quotation  is  taken  from  O.G. Sutton,  Micro-meteorology, 
page  138. 


It  is  suggested  that  Schmidt's  conjugate  power  laws  may  be^more 
accurate.  These  laws  assume  tnat  the  wind  speed  varies  as  the  •=■  power 
of  a  dimensionless  altitude  coordinate  as  it  should  in  idealized 
turbulent  flow,  and  the  diffusivities  vary  as  the  y  power  of  the 


1  1 


altitude  coordinate.  In  two  dimensions,  Sutton  (page  281),  gives  a 
simple  closed  form  solution  for  these  conditions  and  a  source  in  the  form 
of  a  singularity  in  the  x-z  plane.  This  idealization  of  the  source  form 
is  a  frequent  artifice,  used  in  several  fields  of  physics,  e.g.,  electro¬ 
dynamics  and  fluid  dynamics,  whereby  the  singular  vicinity  of  the  source, 
which  produces  the  flow,  is  subsequently  excluded  along  a  streamline, 
equipotential  line,  curve  of  constant  contamination,  or  any  appropriate 
straight  line,  and  only  the  exterior  field  is  considered  as  real. 
Furthermore,  it  is  apparent  that,  if  the  governing  equations  are  linear, 
sums  or  integrals  over  such  sources,  can  be  added  up,  to  give  a  composite 
field,  with  the  desired  properties  (boundaries). 

This  writer  proposes  that  this  “method  of  superposition  of  sources, 
etc."  be  used  to  obtain  the  desired  distribution  of  any  given  contamina¬ 
tion  field.  'ne  may  note  that  the  sources  of  incompressible,  flow  satisfy 
the  Laplace  equation,  which  is  elliptic.  The  "pseudo  sources"  of  super¬ 
sonic  flew  obey  a  hyperbolic  equation  (of  the  type  of  the  wave  equation). 
Since  our  present  sources  obey  the  diffusion  equation,  they  may  be  called 
"parabolic  or  diffusion  sources." 

The  prerequisite  for  such  a  procedure  is  that  several  parabolic 
source  solutions  be  known.  It  is  the  main  objective  in  this  part  to 
derive  several  elementary  three-dimensional  parabolic  source  solutions. 
These  solutions  will  answer  the  question  about  the  distribution  of 
concentration,  when  the  actual  U  and  K  variations  are  the  assumed  ones. 
Beyond  that,  a  superposition  of  such  elementary  solutions  may  enable  us 
in  the  future,  to  deal  with  problems  of  a  more  general  nature. 

THE  CONCENTRATION,  UNDE55.  CONDITIONS  IN  WHICH  U  =  U  •  z7,  Ky  ■  a  •  z7, 

& 

AND  K2  *  b  •  z 7 

The  problem  of  our  three-dimensional  singularity  source  can  be 
rormulated  as  follows: 


We  are  looking  for  a  solution  of  the  parabolic  differential  equation 


0  •  z7  ■  |a 


with  the  following  boundary  conditions: 


K  ^  *  0,  for  z  s  0,  any  x  and  any  y 
q  »  0,  as  x,  y,  or  z— ►  » 


z  =  y  —►0 
1 

U  •  z7  •  q  (x;  y,  z)  =  Q,  independent  of  x 


(2) 

(3) 

(4) 


12 


The  last  condition  represents  the  conservation  of  mass  flow  through  any 
plane  x  -  constant. 


For  a  two-dimensional  singularity  source,  i.e.,  if  in  equation  (1) 
a  is  set  equal  to  zero,  a  solution  has  been  stated  by  Sutton,  page  281, 
which  is  of  the  form: 


(5) 


It  is  readily  verified  by  differentiation,  that  this  solution  with 

.  _  49  U 
A  81  b 

satisiSes  differential  equation  (1)  and  boundary  conditions  (2)  and  (3). 
Condition  (4)  can  be  written,  for  two  dimensions  only,  as: 


and  if 


for  any  fixed  x. 


We  see  from  the  above,  that  M2  is  independent  of  v  and  x,  i.e.,  the 
diffusing  mass  flow  is  constant. 

Sutton's  above  two-dimensional  diffusion  source  is  capable  of 
generalization  to  three  dimensions,  in  several  ways.  To  prove  this, 
let  us  now  consider  the  complete  differential  equation  (1)  and  complete 
conservation  equation  (4). 

Generalizing  "Ansatz"  (5),  set 


whereby  the  coef f icients  p,  m,  and  n  will  be  so  determined,  that  they 
satisfy  the  differential  equation  and  boundary  conditions. 

Inserting  equation  (?)  into  equation  (1),  we  obtain: 
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Next  we  have  to  determine  the  exponents  in  equations  (8),  (9), 

10) ,80  that  they  match.  After  this,  we  shall  match  the  coefficients  of 

;qual  powers  in  z,  y,  and  x. 

'quation  (8)  contains  terms  with  the  following  powers: 

j  /xp+l.  j*“  /xP+2.  Jf  .  y»  /xP« 

Equation  (9)  contains  terms  with  the  following  powers: 

yn-2  .  z7  /XP+1 ;  J  .  y2"-2  /XP+2 

Equation  (10)  contains  terms  with  the  following  powers: 

6  6 

2m-2+7  /xP+l;  2 2m- 2+7  /xP+2 

We  see  that  all  exponents  match,  if  we  set  n  -  2;  and  m  =  y,  as  in  two 
d imens ions . 
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Then: 


Equation  (8)  contains  the  powers 


1  10 
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Equation  (9)  contains  the  powers 

I  1 

z7  /xP+l;  y2  z7  /xP+2 

Equation  (10)  contains  the  powers 

1  10 

z7  /'xP+1 ;  z7  /xP+2 


As  we  see,  p  remains  open. 

Since  we  have  now  matched  all  the  exponents,  we  compare  the 
coefficients  of  equal  powers,  which  provides  th°  following  four  equations 


1 


To  this  add  the  information  that  in  two  dimensions,  if  a  -  B  =  0,  we  know 

that  p  =  Po  *  — ,  so  that  we  have 
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Now  subtracting  equation  (14)  from  equation  (11)  and  setting  p~  +  p^, 
we  get  ^ 


-  P1  U  *  -  a  B  •  2 
U  =  a  B  •  4 

1  25 

which  yields  and  p  «*  p0  +  p^  *  jg 


(11)  -  (14) 
(12) 

(15) 


It  is  noteworthy  that  equation  (15)  holds  true  for  any  values  U,  a,  b. 
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Actually,  these  values  are  determined  by  the  altitude  variations_of 
wind  and  diffusivity.  In  lieu  of  observational  data,  we  set  a  *  b  *  U, 


and  obtain: 

8 

9 

a  72 
*  49 

(14a) 

1 

*  4B 

(12a) 

1 

81  A 

£3  — —  •  A 

49 

(13a) 

Equations  (13a)  and  (14a)  give  the  same  A  value,  and  equation  (12a)  gives 


Therefore,  the  three-dimensional  distribution  of  contamination,  for 
the  indicated  variations  of  U  and  K.,  which  are  in  agreement  with  the 
physical  laws  (of  turbulence)  and  with  observational  data,  assumes  the 
following  simple  form: 


q  (x;  y. 


*18 


exp 


9 

49  z?  i 

81  x  4  x  - 


(16) 


This  is  the  contamination  issuing  from  one  elementary  "parabolic  source," 
mentioned  earlier. 


It  is  obvious  that  solution  (16)  fulfills  not  only  differential 
equation  (1),  but  also  conditions  (2)  and  (3).  But  we  still  have  to 
show  that  it  fulfills  the  condition  of  conservation  of  mass  flow  for 
three  dimensions  (4).  We  can  write  (4),  after  separating  p  into  • 

Po  +  ?L: 


a  Q 


(17) 


Not?,  it  has  been  shown  before,  that  the  first  factor  is  independent  of 
x  and  constant;  it  remains  to  show  that  the  second  factor  is  also 
independent  of  x. 
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S4SS^KC'^55«vr 


-■CEIir  v-  .wi 


We  see  that 


M3  »  M2 


/'  6) 


exp  I  -  B 


lx-?/  J 


Setting  *  r 


M3  “  M2  ‘  /  dr  •  e  _Br*’  =  M2  *  **  •  q.e.d. 


The  diffusing  mass  flow  M3  is  the  same  at  any  plane  x  =  c,  also  in  three 
dimensions . 

6 

THE  SPREAD  OF  CONTAMINATION  q  (x;  y,z),  FOR  U  =  U  •  z7, 

6  6 
Ky  =  a  •  z^,  Kz  =  b  •  z7 

Another  chree-dimensional  singularity  source,  which  has  been 
suggested  as  representing  actual  meteorological  conditions,  is  the  one 
indicated  above,  so  that  we  have: 

“  -  *7 '  12 "  §  (a  •  *7 '  I?)  +  it  (b  •  J  ■  It) 


Setting  as  before: 


q  =  £_  .exp 


we  get  this  time: 
6 


-  A  •  2^  -  B  •  exp  [-  L  (x,y,z)| 

x  x  xP  1  J 


l-  U  ‘  zT  •  f?  =  u  ’  zl  '  +  +  *  exP  [-  L  (x»y»z)] 

6  0  6  o  3 

o  7  dzq  7  B*C  ,  .>  n_2  ,  B^‘ C  2  2n-2 1 

2.  a  •  z'  •  =  a  •  z!  -  n(n-l)  ♦  yn  *■  +  n^  •  y^u  *•  j  • 

•  exp  j-L  (x,y,z)j 

3.  J.L  .  J  .  gal  .  f.  bA-.CjgL.l_t  |)  z(">-2+f)+  bA^m2  .  z(2m-2+4 


|-(b  J .  fa 

dz  dz 


exp  J-L(x,y,z)J 
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1  6.i_.  6 

1.  Contains  the  powers  z?/xP+~;  zJ  /xP^2;  z7  .  ,.n/xp+2 

-  6 

2.  Contains  the  powers  z?/xP+*;  z?  •  y2n-2/xp+2 

_6  *> 

3.  Contains  the  powers  z(m~2)+7/xp+l.  2(2m-2,»+7/xp+2 

All  powers  will  be  matched,  if  we  take  n  <*  2,  m  =  2,  and  leave  p  open, 
at  first. 

By  comparing  the  coefficients  of  equal  powers  we  get  the  equations: 

J& 


7/_ _ 

xP+1 

|+2 

z' 

^FT 

6  2 
zi  •  y2 


-  U  •  p  =  -  a  *  B  •  2  -  b  *  A  ■  2 


U  =  b  A  *  4 


13 


U  -  a  B  •  4 


(X) 

(ID 

(III) 


In  two  dimensions,  where  a  =  B  =  0,  we  get 


U  ‘  Po  “  b  A  y 

U  ■  b  A  •  4 
Hence: 

26  1  ^  I 

4p0  *  ~j\  po  =  for  any  U  and  b  value. 

Set  again,  p  «  pQ  and  p^,  then 
-  U  •'  P|  =  -  a  B  •  2 
U  =  a  B  •  4 

4px  *  2,  px  «  j  again. 

If  now  we  set  U  =  b  *=  a,  then 


(IV) 

(in 


(i)  -  (tv) 
(in) 


THE  CASE:  U  -  B 


K- 


K, 


b  •  z 


u 


2. 


Ck-, 

dx 


32q 


3  z  ' 


^pTT  +  Xp-KT  4  *  exP  ["  L  (x>y»g)j 

f  B* S  ,  n-2  ,  B2-C  ?  2n-?l 

2  {“  ~~fY  n(n-l)  *  yn  2  +  -^p?  r.4  •  y*n  2j  • 

l  x‘  xt* 

f  .1 

•  exp  -  «,  U»y*s)| 


> 

3. 

i|‘- 

,  ia 

*z 

1. 

Contains 

the 

) 

2. 

Contains 

the 

3 . 

Contains 

the 

bA*C*m^ 


“x^ 


,tn-i 


bA2*C‘m2 

*  2  j 


exp  -  L  (x,y,z 


We  can  match  the  powers,  if  we  set  m  =  2  and  n  **  2,  for  three  dimensions. 
For  the  special,  two-dimensional  case  we  set  a  =  B  1  0  and  n  =»  0.  Then 


-U  •  p0  -  -  b  A  •  4 
U  »  b  A  •  4 

=  2.  For  three  dimensions 
b  A  ’  4;  U  ■=  a  B  *  4 
Hence  4p-,  »  2,  and  =  i.,  again. 


Therefore  pQ  »  1  with  n  -  0,  m 
-U  •  p1  =  -  a  B  •  2;  U  = 


The  addition  of  the  third  dimension,  appears  always  to  increase  the 
negative  power  of  x  by  one  half,  so  long  as  the  crosswind  distribution 
of  q,  i.e.,  Ky  is  not  made  to  depend  on  y. 

THE  CASE:  U  (z)  =  U  •  z;  Ky  *  a  z  *  y;  K»  *  b  *  z 

The  meteorological  data  seem  to  require  that  the  fall-off  of 
concentration  in  the  downwind  direction  is  somewhere  between  x~*  and 
x~^‘5.  Since  the  preceding  cases  seem  to  indicate  that  any  z-dependence 
of  Ky  does  not  lead  beyond  p  *  1.5,  we  try  to  gain  new  insight,  by  making 
Ky  “  a  •  z  ’  y. 
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While  equations  (1)  and  (3)  of  case  5  remain  the  same,  we  have  now 
to  replace  equation  (2)  by 


Therefore  the  terms  (1),  (2a),  and  (3)  of  our  differential  equation 
contain  the  following  powers  of  x,y,  and  z: 

1 .  z  / xP+*  ;  z13"*^  / xP+2 ;  z  •  y R/ xP+- 

2a.  z  •  yn~’/xP+-*-;  z  •  y2n-l/xp+2 

3.  zflJ-l/xP+i;  z2m-l/xp+2 

We  see  that  we  can  match  all  powers  by  setting  n  »  1,  m  =  2.  First, 
solve  the  two-dimensional  case,  by  setting  a  =  B  =  0,  m  =  2.  Then, 
as  in  case  3 

-  U  '  pQ  =  “  b  A  •  4;  'u  =  b  A  •  4.  Hence  pG  =  l,  again. 

For  three  dimensions: 

•  U  •  p  =  •  a  B  *  1  -  b  A  •  4;  U  s  b  A  ■  4;  U  *  a  B  ■  1 

Hence  p  -  pQ  =  =  -f  1  and  p  =  2.  The  fall-off  of  q(x,0)  is  propor¬ 

tional  to  x"-. 
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i_ 

THE  CASE:  U(z)  =  U*z;Ky  =  a-  z-y2;K2«b-z 

o’ 

Since  the  fall-off  in  case  6  was  too  large,  we  try  Ky  *  a  •  z  •  yA. 
Again  the  terms  (1)  and  (3)  remain  the  same  as  in  cases  5  and  6.  But 
term  (2)  must  be  replaced  by  the  following: 


d_ 

dy 


azy 


,2  23. 


oy 


=  j  a  2  •  y 


i 


dq  Y  d^q 

dy  dy 


(2b) 


Substituting  expression  (7;  for  q,  we  get: 


a  z  *  y 


.  v  2  .  2a  » I 


dy  2 


a  z  *  y 


i. 

'2 


B  C-n-y0-1 


•  exp 


-  L  (x,y,z) 


jL  « 


+  a  z  yA  *  -azy 


1 

,2 


B^C« R2 .y2n-2  ^  BC'n(n-l) -yn~2  | 


X' 


pTT 


exp  -L 


Hence: 


3 

3 

B*C*n-yn“2 

\  +  n-1 

■ C'n^*y^n"2 

•  exp 

-  L  (x,y,z) 

“  xP+I 

+  —^prL 

(2b) 


The  terms  (1),  (2b),  and  (3)  of  our  differential  equation  contain  the 
following  powers  of  x,y,  and  z: 


1.  z/xP+1-;  z1^1  / xP'i'2  j  z  .  yn/xp+2 

2b.  z  •  yD~^/xP+-*- ;  z  •  y2n-^yxp+2 
3.  zm-l/xp4'l.  s2m-l/xp+2 

O 

All  terms  can  be  matched,  if  we  take  m  =  2,  n  =  From  the  two- 
dimensional  case,  (a  =  B  *»  0)  ,  we  get  as  before: 


Po  =  1 

For  three  dimensions,  we  get  the  equations: 

-U-p*-aB'-|-bA-  4;  U  =  b  A  «  4;  U  =  a  B  ■  | 

Subtracting  from  the  first  equation;  -U  *  pQ  =  -  b  A- •  4,  we  obtain  now 
P  “  P0  =  Pi  =  f »  80  that  P  =  1-66 

This  fall-off  of  q(x,0)  seems  to  agree  fairly  well  with  the  observational 
data. 


If  necessary,  new  parabolic  source  solutions  can  be  readily 
generated . 
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CONCLUSIONS 


We  have  shown  that,  depending  on  the  variation  of  the  diffusivities 
Ky  and  K  with  the  altitude  coordinate  t  and  crosswind  coordinate  y,  the 
rate  of  decline  of  concentration  q,  with  downwind  distance,  xP  may 
follow  any  power  law  between  p  ~  1.5  and  p~  2.0.  In  the  crosswind 
direction,  a  Gaussian  distribution  results.  This  general  behavior  of 
the  distribution  of  concentration,  at  ground  level,  is  well  substantiated 
by  observations  in  all  kinds  of  flat  lands,  as  reported  in  the  literature. 
An  empirical  exponent  which  is  favored  is  p  —  i -75*  * 


*  Pasquil,  F.  Atmospheric  Diffusion,  London,  Van  Nostrand,  Ltd., 
1962.  Page  83. 


